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ABSTRACT 

Fluctuations in the brightness of the background radiation can lead to confusion with 
real point sources. Such background emission confusion will be important for infrared 
observations with relatively large beam sizes since the amount of fluctuation tends to 
increase with angular scale. In order to quantitively assess the effect of the background 
emission on the detection of point sources for current and future far-infrared obser- 
vations by space-borne missions such as Spitzer, ASTRO-F, Herschel and SPICA, we 
have extended the Galactic emission map to higher angular resolution than the cur- 
rently available data. Using this high resolution map, we estimate the sky confusion 
noise due to the emission from interstellar dust clouds or cirrus, based on fluctuation 
analysis and detailed photometry over realistically simulated images. We find that 
the confusion noise derived by simple fluctuation analysis agrees well with the result 
from realistic simulations. Although the sky confusion noise becomes dominant in long 
wavelength bands (> 100 /im) with 60 - 90cm aperture missions, it is expected to be 
two order of magnitude smaller for the next generation space missions with larger 
aperture sizes such as Herschel and SPICA. 

Key vifords: methods: data analysis - techniques: image processing - ISM: structure 
- galaxies: photometry - Infrared: ISM 



1 INTRODUCTION 

The detection of faint sources in far IR can be greatly af- 
fected by the amount and structure of the background ra- 
diation. The main source of background radiation in far IR 
is the smooth component of the Galactic emission, known 
as cirrus emission. The amount of emission manifests itself 
as photon noise whose fluctuations follow Poisson statistics. 
In addition, any brightness fluctuation at scales below the 
beam size could cause confusion with real point sources. The 
cirrus emission was discovered by the Infrared Astronomy 
Satellite {IRAS) JLow et ai. 1984(l . and is thought to be due 
to radiatively heated interstellar dust in irregular clouds of 
wide ranges of spatial scales. The cirrus emission peaks at 
far-IR wavelengths but was detected in all four IRAS bands 
at 12, 25, 60, and 100 ^m (Helou & Beichman 1990, hereafter 
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HBO). The brightness of cirrus emission depends upon the 
Galactic latitude and is significant for wavelengths longer 
than 60 /im. The cirrus emission, which is the main source of 
background radiation in far-IR, causes an uncertainty in the 
determination of source fluxes use its brightness varies from 
place to place. The accurate determination of observational 
detection limits requires a knowledge of the cirrus emission 
as a function of position on the sky. The other important 
factor affecting the source detection is the source confusion 
which mainly depends upon the telescope beam size and the 
source distribution itself. The effects resulting from a com- 
bination of the sky confusion and the source confusion will 
be discussed in depth in the forthcoming paper [Jeong et al. 
2004c {Paper lT| l, in preparation], and we concentrate on the 
effect of sky confusion in the present paper. 

There have been realistic estimations of the sky confu- 
sion from observational data from IRAS and the Infrared 
Space Observatory {ISO) (Gautier et al. 1992; HB90; Herb- 
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Figure 1. Schematic outline of the reference aperture configura- 
tions for two symmetrically placed circular apertures (Gautier et 
al. 1992). 



stmeier et al. 1998; Kiss et al. 2001). However, the reso- 
lution of the data from IRAS and ISO is not sufficient 
to the apphcation to larger missions planned in future. 
Many valuable data in the far-IR wavelength range will be 
available within or around this decade by a multitude of 
IR space projects such as S pitzer ([Gallagher_et_aL^^ 
ASTRO-F (|M urakami, T998l IShibai 20001 |Nakagawa^OOl| 
IPearson et al. 2004). Herschel Space Observatory (HSO) 
IjPilbratt 2003 Poglit sch et al. 2003^ and the Space In- 
frared Telescope for Cosmology and Astrophysics (SPICA) 
| |Nakagawa 2004| l. Since these instruments will observe the 
sky with high sensitivities and high angular resolution, it is 
necessary to understand the factors determining their detec- 
tion limits. 

The purpose of the present paper is to investigate the ef- 
fects of cirrus emission on the detection of faint point sources 
in highly sensitive future infrared observations. Based on the 
measured power spectrum and the spectral energy distribu- 
tion models of the dust emission over the entire sky, we 
generate the dust map with higher spatial resolution in var- 
ious relevant wavelength bands by extrapolating the power 
spectrum to small scales. 

This paper is organized as follows. In Section |5| we 
briefly describe the sky confusion noise due to sky bright- 
ness fluctuations. In Sectional the high angular resolution 
realization of Galactic dust emission in various IR bands is 
presented. Based upon the speciflcations of each IR mission, 
we estimate the sky confusion noise by using simple fluctua- 
tion analysis in Section |1] We compare estimated detection 
limits based on fluctuation analysis with the results based 
on the photometry on realistically simulated data in Section 
13 Our conclusions are summarised in Section |S| 



2 CONFUSION DUE TO SKY FLUCTUATION 

Measuring the brightness of sources involves subtracting the 
sky background derived from the well-deflned reference. The 
fluctuations in the surface brightness of extended structure 
on similar scales to the resolution of the telescope and in- 
strument beam can produce spurious events that can be eas- 
ily mistaken for genuine point sources. This is because the 
source detection is usually simply accomplished from the dif- 
ference in signal between the on-source position and some 
background position. Therefore sky confusion noise due to 
the sky brightness fluctuations, N{9), is defined as (HB90; 
Gautier et al. 1992): 



where f2 is the solid angle of the measuring aperture, 6 is 
the angular separation between the target and reference sky 
positions, and S{6) is the second order structure function, 
which is defined as tGautier et al. 1992t : 



S{e) = 



I(a 



i{x-e) + i{x + ( 



(2) 



where / is the sky brightness, x is the location of the target, 
and ( ) represents the average taken over the whole map. 
For the configuration of two symmetrically placed reference 
apertures, see Fig. Q 

Although the zodiacal emission is main background 
source in the short wavelength of far-IR range in low eclip- 
tic latitude regions, it will not contribute to the fluctuations 
on the large scales because the zodiacal fight is generally 
smooth on scales smaller than typical resolution of IR ob- 
servations IIReach et al. i ggS'l^ safi et al. 199811 . From the 
analysis of the ISO data, Abraham et al. 1199711 searched 
for the brightness fiuctuations in the zodiacal light at 25 
/im with 5 fields of ~ 0.5° x 0.5° at low, intermediate, and 
high ecliptic latitudes. They found that an upper limit to 
the fiuctuations of 0.2 per cent of the total brightness level 
was estimated for an aperture of 3' diameter. This amount 
of fiuctuations would not cause any significant noise. 

Therefore, the sky confusion noise is mainly related to 
the spatial properties of the cirrus. In many cases, the power 
spectrum of the dust emission can be expressed as a simple 
power-law. Using the IRAS data at 100 pm, Gautier et al. 
I|1992|l computed the power spectrum P of the spatial fluc- 
tuations of cirrus emission as a function of spatial frequency 
k, for angles between 4' and 400'. 



(3) 



where d represents the angular scale corresponding angular 
frequency {k — 2ZL), The subscript on and d denotes a 
reference scale, Pq is the powers at k = ko, and a is the 
index of the power spectrum. Since the second order struc- 
ture function is proportional to power spectrum representing 
the spatial structure of cirrus, the sky confusion noise A'^ on 
a scale d corresponding to the width of the measurement 
aperture scales as: 



N cx 



(do) 



(4) 



HB90 extended the work by Gautier et al. I|1992^ at A = 
100 fim in order to estimate the sky confusion at all wave- 
lengths, using the empirical relationship, Po oc (7o>^ and 
a = —3 in Gautier et al. (119921 . They found an approxima- 
tion for the cirrus confusion noise as follows (hereafter HB90 
formula) : 



A 



100 fim 



Dt 
1 m 



(A) 



1 MJysr-i 



mjy, (5) 



N{e) = n./s{e). 



(1) 



where is a constant, A the wavelength of the measure- 
ment, Dt the diameter of the telescope, and (/a) is the mean 
brightness at the observation wavelength. They set the con- 
stant C, to be 0.3. 

This indicates that the sky confusion depends upon 
both the variation of the surface brightness in the back- 
ground structure and the resolution of the telescope. Conse- 
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Figure 2. Measured power spectrum of dust emission in the dust 
map by SFD98 (Schlegel, Finkbeiner & Davis 1998). The four 
curves represent four patches selected in the Northern and the 
Southern Galactic sky at fe = |50|°. 



quently, the noise becomes less significant for larger aperture 
sizes. 



3 GENERATION OF CIRRUS MAP 

In order to investigate the sky confusion for the present and 
upcoming infrared space missions with a high resolution, we 
need the information on the behavior of cirrus emission in 
very small scales. Since observationally available data have 
rather low resolution, we need to add high resolution com- 
ponent. In this section, we describe the method of extending 
the low resolution data to high resolution. For the observa- 
tional low resolution data, we used the all-sky 100 /im dust 
map generated from the IRAS and COBE data by Schlegel, 
Finkbeiner, and Davis (1998; hereafter SFD98). 



3.1 Fluctuations at Higher Spatial Resolution 

3.1.1 Measured Power Spectrum 

Fig.|5|shows the measured power spectrum in the dust maps 
of SFD98 at a Galactic latitude of & = |50| degrees. These 
power spectra are well fitted to power laws of index -2.9. 
However, the power drops at higher frequencies correspond- 
ing to the map resolution of ~ 6.1 arcmin. This breakdown 
of the power spectrum is due to the large beam size of IRAS 
map. Although we can recover the small-scale fluctuation by 
the deconvolution of a point spread function (PSF), there is 
clearly some limitation. We need to generate the dust map 
including the contributions from small-scale fluctuations in 
order to study for the planned present and future missions 
with high resolution (< 1 arcmin). We obtain such high reso- 
lution map by adding small-scale structure of cirrus emission 
to the low-resolution map of SFD98 assuming that the small- 
scale fluctuations also follow the estimated power spectrum 
with the same power-law index, as described above. 



3.1.2 Small Scale of Fluctuations 

The power, P{k), is defined as the variance of the amplitude 
in the fluctuations: 

P(k) ^ (I S, f) = ^J iixf^^Anx'dx, (6) 

where 5k is the perturbation field, (| 5k P) is the variance 
of the fiuctuation and ^(x) is the correlation function of the 
brightness field. We assume that the distribution of fluctua- 
tions is approximated as a random Gaussian process where 
the Fourier components 5^ have random phases so that the 
statistical properties of distribution are fully described by 
the power spectrum | 5k P (IPeebles 198011 . In this case, 
we can set each fluctuation within a finite grid in the fre- 
quency domain by a random Gaussian process of the am- 
plitude of each fluctuation considering the realization of a 
volume for the sample embedded within a larger flnite vol- 
ume lIGott et al. 19901 IFar^et al. 1994| IPeacock 1999t . We 
assign Fourier amplitudes randomly within the above dis- 
tribution in the flnite volume and assign phases randomly 
between and 2-n. Since the field used in this simulation is 
small (< 10 degree), we can take the small-angle approxima- 
tion and treat the patch of sky as flat I White et al. 19991 . In 
the flat sky approximation, we obtain the power spectrum 
and generate a patch of the dust map in cartesian coordi- 
nates. 

We generate a realistic distribution of the Galactic emis- 
sion in the following manner. The basic data for the infor- 
mation of the large-scale structure are obtained from the 
low resolution all-sky map by SFD98. We add the simulated 
small-scale structure to these basic data in the Fourier do- 
main, where the power spectrum of the small-scale structure 
follows that of the large-scale structure. Fig. |3 shows our 
simulated emission map including small-scale fluctuations. 
The left panel of Fig. El shows the simulated dust emission 
image corresponding to a power spectrum with a = —3. The 
middle panel includes only the emission above the resolution 
of the dust map by SFD98, ~ 6.1 arcmin, (large-scale emis- 
sion) while the right panel shows the emission above the 
resolution of the dust map by SFD98 (separated in Fourier 
domain, i.e., small-scale emission). The lower panel shows 
the profiles for selected areas of two images (upper-left and 
upper-middle panels). We find in this simulation that the 
emission including the high resolution, small-scale compo- 
nent (above the resolution of the dust map by SFD98 to a 
resolution of 4 arcsec) reflects the trend of the large-scale 
emission (above the resolution of SFD98 dust map). 

We obtain a patch of the dust map including small- 
scale fluctuations by summing the large-scale component of 
SFD98 map and the small-scale component of the simulated 
emission in the Fourier domain. According to this scheme 
of Fourier power spectrum analysis, the cutoff spatial fre- 
quency of the dust map by SFD98 is set to the Nyquist limit, 
i.e. a half of the spatial frequency corresponding to the res- 
olution of the dust map by SFD98. We use the power spec- 
trum fitted below the Nyquist sampling limit in order to ex- 
tend the power spectrum to higher spatial frequencies. Typ- 
ically, the 2D power spectrum of a SFD98 dust map patch 
shows the presence of a cross along spatial frequencies of x 
and y axis if we assume that the centre in the spatial domain 
is regarded as the spatial frequency 0. This cross is caused 
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Figure 3. Simulated dust emission map (upper) and the profile of map (lower). The upper-left panel shows the simulated image assuming 
a power spectrum with a power index of -3. The upper-middle panel and the upper-right panel show only large-scale fluctuations and 
small-scale fluctuations, respectively. The lower panel shows the one-dimensional profile for a selected part of the upper-left and the 
upper-middle panel. 



by the Fast Fourier Transform (FFT) algorithm that makes 
an "infinite pavement" with the image prior to computing 
the Fourier transform IMiville-Deschenes et al. 2002t . In or- 
der to preserve the information of the emission at the edges, 
we directly use the power at the spatial frequencies of x and 
y axis, and extrapolate the power at other spatial frequen- 
cies (above the cutoff spatial frequency) according to the 
estimated power spectrum. In Fig. 0] we show a patch of 
the dust map by SFD98 at a Galactic latitude of 50 degree 
(upper left), a patch regenerated by extending the power 
spectrum (upper right) and the estimated power spectrum 
(lower panel). 



3.2 Dust Emission at Other Wavelengths 

Assuming that the spatial structure of the dust emission is 
independent of wavelength, we can obtain the dust map at 
other wavelengths than 100 ^m by applying an appropriate 
model for the Spectral Energy Distribution (SED). Since 
the dust particles are small (< 0.25 fim) compared with 
far-IR wavelengths, the opacity does not depend upon the 
details of the particle size distribution, but on the nature 



of the emitting material itself. In the far-IR, the opacity 
generally follows a power law: 

K„ (X 1/'^ (7) 

with frequency v. 

The SED may be approximated as one-component or 
two-component models jSchlegei, I^'inkbeiner fc Davis 1998| 
IFinkbeiner et al. 199911 . The dust temperature map is 
constructed from the COBE Diffuse Infrared Back- 
ground Experiment (DIRBE) 100 /im and 240 /^m data 
l |Boggess et al. 1992| l which was designed to search for the 
cosmic IR background radiation. For a one-component 
moedel, the emission I^, at frequency i/ can be expressed 
as 



-^ioo(/5i^)-'^ioo 



(8) 



where Bv{T) is the Planck function at temperature T, Jioo is 
the Z)/i?_BiJ-calibrated 100 /im map, -R'j^q (/3, T) is the colour 
correction factor for the DIRBE 100 filter when observ- 
ing a Bv{T) spectrum {DIRBE Explanatory Supplement 
1995). Although the generated temperature maps have rel- 
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Figure 4. Patch of SFD98 dust map, regenerated patch (upper panel) and the estimated power spectrum (lower panel). The upper-left 
panel is a patch of the SFD98 dust map at the Galactic latitude of 50 degree and the upper-right panel is the regenerated patch based 
upon the patch from the SFD98 dust map. The dashed and solid lines in the lower panel show the estimated power spectrum of the 
upper-left and the upper-right panels, respectively. Note that the Nyquist frequency in the power spectrum of the upper-right panel is 7.5 
arcmin"^, but we only plot to ~ 0.5 arcmin"^. The dotted line shows the fit to the power spectrum below the spatial cutoff frequency. 



atively low resolution (1.3°) compared with our simulated 
dust map patch, we interpolate this map to small grid sizes 
(< 10 arcsec). Taking the emissivity model with /3 = 2 
JDraine fc Lee 198411 . we can obtain the dust temperature 
from the DIRBE 100 /im/240 /im emission ratio. 

Based upon laboratory measurements, a multicompo- 
nent model for interstellar dust has been constructed by 
Pollack et al. 119941 . In order to solve the inconsistency 
of the u'^ emissivity model in the 100 - 2100 GHz (3000 
— 143 ^m) emission, Finkbeiner et al. 1199911 used a two- 
component model where diverse grain species dominate the 
emission at different frequencies in order to fit the data of the 
COBE Far Infrared Absolute Spectrophotometer (FIRAS). 
Assuming that each component of the dust has a power-law 
emissivity over the FIRAS range, Finkbeiner et al. (I1999II 
constructed the emission 1^ in multicomponent model: 



h = 



(9) 



where fi is a normalization factor for the i-th grain com- 
ponent, Ti is the temperature of component i, Kioo is the 
DIRBE colour-correction factor and /loo is the SFD98 100 
^m flux in the DIRBE filter. The emission efficiency Qi{v) 
is the ratio of the emission cross section to the geometrical 
cross section of the grain component i. In order to obtain 
the temperature of each component, we further need effec- 
tive absorption opacity deflned by 



(10) 



where k^''^ is the absorption opacity of the i-th component, 
and JisRF is the mean intensity of interstellar radiation field. 
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Figure 5. Comparison between the one-component dust model 
and tile two-component dust model for one small patch. The dust 
emission of the two-component model in the wavelength range 
from 120 /^m to 200 ^m is slightly higher than that of the one- 
component model due to the dominant contribution by carbon 
grains. 



Finkbeiner et al. Ill999t assumed that the normalization fac- 
tors do not vary with locations and size independent optical 
properties of dust grains. The emission efficiency factor Qi 
at far-IR is further assumed to follow a power-law with dif- 
ferent indices (/3) for different dust species. In the present 
work, we adopted the 'best-fitting' two-component model by 
Finkbeiner et al. (19991: /3i = 1.67, 132=2.70, fi = 0.0363, 
/2 = 0.9637, and (71/52 = 13.0, where qi = Kf%uo)/K* 
which represents the ratio of far-IR emission cross section 
to the UV/optical absorption cross section. The reference 
frequency uq is that corresponding to wavelength 100 /im. 

If we further assume that the interstellar radiation field 
has constant spectrum, the temperature of each compo- 
nent can be uniquely determined by the far-IR spectrum 
represented by the DIRBE 100 /im/240 /im ratio. A two- 
component model provides a fit to an accuracy of ~ 15 per 
cent to all the FIRAS data over the entire high-latitude sky. 
In Fig. 13 we see the dust emission for the one-component 
and two-component dust models [see Schlegel et al. 1199811 : 
Finkbeiner et al. (19991]. The two-component model agrees 
better with the FIRAS data in the wavelength range longer 
than 100 ^ra where the dust emission estimated from one- 
component model is significantly lower than the estimate 
from the two-component model. 

In two models, the contribution of the small grains re- 
sulting in an excess below 100 /im is not considered. Since 
there is no significant difference between models below 100 
(Um while the dust emission of the two-component model 
is slightly higher than that of the one-component model in 
wavelengths ranging from 120 to 200 /im, we use the two- 
component model in our calculations. 

Through a PSF convolution at each wavelength and a 
wavelength integration over a 5 /im wavelength grid, we ob- 
tain the high resolution dust map in other bands. 



Table 2. Sky confusion noise estimated from HB90 formula for 
each space mission. The instrumental parameters for each mission 
are given in Table [Tl The mean brightness here is fixed to be 1 
MJy sr-i. 



N (mjy) 



Space Mission 


SW 


LW 


ISO 


0.83 


4.05 


Spitzer 


0.18 


1.46 


ASTRO-F 


0.40 


1.89 


Herschel 


0.0054 


0.042 


SPICA 


0.0054 


0.042 



4 FLUCTUATION ANALYSIS FOR SKY 
CONFUSION NOISE 

Among the parameters affecting the sky confusion noise, 
most of them depend upon the mean brightness, the spatial 
structure of the cirrus, and the observing wavelength, as seen 
in equation l^. In TableQ we list the basic instrumental pa- 
rameters of present and future IR space missions; the aper- 
ture of the telescope, Full Width at Half Maximum (FWHM) 
of the beam profile and the pixel size for each detector. For 
comparison with previous studies (iHerbstmeier et al. 19981 
IKiss et al. 200111 . we include the specifications for ISO. We 
select a short wavelength band (SW) and a long wavelength 
band (LW) for each mission. 

In order to examine the dependency of the sky confusion 
noise on the instrumental parameters, we list sky confusion 
A*' estimated from HB90 formula for each mission considered 
in this work in Table |5| As the aperture of the telescope 
becomes larger or the wavelength becomes shorter, sky con- 
fusion A'^ should become correspondingly smaller. In Section 
13 we obtained the dust maps extended to high spatial reso- 
lution over a wide spectral range. With this simulated dust 
map, we estimate the sky confusion noise for various space 
mission projects. 

4.1 Selected Regions 

We generate the PSF-convolved patches of a dust map as 
a function of increasing Galactic latitude (decreasing sky 
brightness) from 0.3 MJy sr~^ to 25 MJy sr~^ at 100 fim 
at a resolution of 1 arcsec by using the method explained in 
Section 13 The size of the simulated image is 1.3° x 1.3°. 
For the PSF, we used an ideal circular aperture Airy pat- 
tern corresponding to the aperture size of telescopes. In Fig. 
13 we can see the PSF-convolved small patch of dust map 
(900" X 900") for each space mission. As the aperture of 
the telescope becomes larger, the structures that can be vis- 
ible become smaller. Since the cirrus emission generally de- 
pends upon the Galactic latitude, we select the patches as 
a function of the Galactic latitude. We list the properties 
of selected regions at a Galactic longitude of 0° among 50 
patches in Table |3 The estimated power spectrum in Table 
13 differs from patch to patch. In order to reflect the large 
structure of the dust map and reduce the discrepancies of 
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Table 1. Instrumental parameters for various space missions. 



Aperture Wavelength FWHM " Pixel size 
(meter) (m™) (arcsec) (arcsec) 



Space Mission 




SW 


LW 


SW 


LW 


SW 


LW 


ISO 


0.6 


90 


170 


31.8 


60 


46 


92 


Spitzer 


0.85 


70 


160 


16.7 


35.2 


9.84 


16 


ASTRO-F 


0.67 


75 


140 


23 


44 


26.8 


44.2 


Herschel 


3.5 


70 


160 


4.3 


9.7 


3.2 


6.4 


SPICA 


3.5 


70 


160 


4.3 


9.7 


1.8 


3.6 



" FWHM of diffraction pattern. 

>> Two ISOPHOT filters (C1.90 in SW band and C2.170 in LW band). 
MIPS bands for the Spitzer mission. 

ASTRO-F/FIS (Far Infrared Surveyor) has a WIDE-S band in SW and WIDE-L band in LW. 
' PACS have 'blue' array in short wavelength (60-85/im or 85-130^m) and the 'red' array in long wavelength (130-210^m). 




Figure 6. PSF-convolved patch of the dust map for space mis- 
sion; ISO (upper-left), ASTRO-F (upper-right), Spitzer (lower- 
left), Herschel/SPICA (lower-right) missions. 

the power spectrum between adjacent patches, we use a large 
area around the patch (~ 2.5° x 2.5°) in the measurement 
of the power spectrum. 



Table 3. Properties of the selected regions. The Galactic lon- 
gitude of all patches is 0°. Iq is a mean sky brightness, a is 
the power index of the power spectrum, and Pq is the power 
estimated at 0.01 arcmin"^ and 100 /.tm. 







lo " 






log Po 






(MJy sr" 






(Jy2 sr-l) 


Region " 


70/^m 


lOO^tm 


160/im 






b=10° 


5.4 


24.4 


53.9 


-3.45±0.11 


9.00±0.17 


6=17° 


3.5 


18.6 


45.3 


-3.50±0.16 


9.05±0.24 


6=22° 


3.5 


15.3 


34.1 


-3.54±0.15 


8.48±0.22 


6=28° 


2.2 


8.9 


24.7 


-3.50±0.15 


7.74±0.21 


6=36° 


1.2 


6.0 


14.4 


-3.80±0.10 


7.41±0.15 


6=45° 


0.6 


2.8 


6.2 


-3.13±0.12 


6.39±0.18 


6=59° 


0.3 


1.4 


2.9 


-2.99±0.09 


6.00±0.13 


6=70° 


0.2 


1.2 


2.6 


-3.20±0.10 


6.27±0.15 


6=84° 


0.1 


0.8 


1.8 


-2.87±0.09 


5.77±0.14 


6=90° 


0.1 


0.5 


1.4 


-2.87±0.08 


5.66±0.12 



4.2 Estimation of Sky Confusion Noise 

4-2.1 Contribution of Instrumental Noise 

In order to estimate the sky confusion noise, the 
structure function for the cirrus emission patch ob- 
tained by measuring the sky brightness fluctuations is 
widely used (jGauticr et al. 1992 Herbstmcicr et al. 1998 
IKiss et al. 200l|l . The size of the measuring aperture is set to 
be the FWHM of each beam profile if the detector pixel size 
is smaller than the FWHM of a beam profile. Since the sky 
confusion noise and the instrumental noise are statistically 



independent ijHerbstmeier et al. 19981 IKiss et al. 2001^ . the 

measured noise A'^mcas is 

^'"mcas = A^^ + ?? • Cr^t, (11) 

where A'' is the sky confusion noise corresponding la, (Tinst 
is the instrumental noise, and rj is the contribution factor 
from the instrumental noise. The contribution factor r] can 
be determined by the size of the measurement aperture and 
the separation (see equation |5| and Fig. 0. 
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Figure 7. The fraction of the sky brightness for all sky. Note 
that most of the sky have the sky brightness below 1 MJy sr~^ 
(SW) and 15 MJy sr'^ (LW). The contribution in the highest 
mean brightness resulted from near the Galactic center. 



4-2.2 Comparison with Other Results 

We estimate the sky confusion noise from the patches of the 
simulated sky map. In Fig.|2l we plot the fractional area as a 
function of sky brightness over the whole sky to visualize the 
sky brightness distribution. Since we consider the sky con- 
fusion caused solely by the emission from cirrus structures, 
we do not include any contribution from the instrumental 
noise. 

In order to determine a dependency of the sky confu- 
sion noise on separation, we performed a "calculation" for 
the estimation of sky confusion noise for given mean bright- 
ness of the sky patch for each space mission {ISO, Spitzer, 
ASTRO-F, and Herschel/SPICA) by systematically vary- 
ing the value of s from 2 to 7, using equation where s 
parameter is related to the separation 9 = sD. Generally, 
larger separation causes larger sky confusion noise because 
we may be estimating the fluctuations from different struc- 
tures. In practical photometry, large separations are gener- 
ally used, i.e., 6 = sD, s > 2 in the configuration of Fig. 
m iKiss et al. 20011 [Laureijs et al. 2003| l. As a reference, we 
take the estimate of the sky confusion noise with s = 2.5 
for a comparison of the measured sky confusion with the 
photometric results given in Section |S] In the source detec- 
tion, the background estimation parameter have the same 
role with the separation parameter. We found the optimal 
value for the background estimation parameter through the 
photometry (see Section fS. 21 for detailed explanation). 

In Figs |H| - 111! we present our estimates of the sky 
confusion noise for the ISO, Spitzer, ASTRO-F and Her- 
schel/SPICA space missions comparing the formula for the 
sky confusion noise predicted by HB90 (hereafter HB90 for- 
mula). For ISO results, the sky confusion noise with s = 2.5 
is overestimated for the dark fields, but underestimated 
for the bright fields (see Fig. |HJ. With larger separations, 
e.g., s = 7, the estimated confusion noise approaches the 
HB90 formula although it is still overestimated for the dark 
fields. We can see the same tendency in other studies in 
the sky confusion noise measured from ISO observations 
IjHerbstmeier et al. 19981 IKiss et al. 200l| l. The measured 
sky confusion noise for the Spitzer and Herschel/SPICA mis- 



sions are much lower than the predictions of HB90 except 
for the dark fields (see Figs 1101 and II H . 

Comparing the empirical relation between Pa and Jo by 
Gautier et al. (1992 ), we present our estimated Po in Fig. 
1121 It shows a lower Po in bright fields and the higher Po 
in dark fields could cause an underestimation in the bright 
fields and an overestimation in the dark fields of the sky 
confusion noise. Such inconsistencies, overestimation of Po 
in bright fields and underestimation of Po in dark fields, also 
appear in other regions of the sky. By fitting our estima- 
tions of Po, we obtained a new relation between the Po and 
7o. The HB90 formula assumed the wavelength dependency 
only through the beam size. However, although the cirrus 
structure is generally preserved in other wavelengths, the 
empirical relation should be scaled according to the vari- 
ation of the cirrus brightness with wavelength, i.e, cirrus 
spectral energy distribution. Therefore, in order to apply 
our empirical formula to other wavelength bands, we need 
some additional correction. For this correction, we used the 
ratio of the mean brightness at the two wavelengths, e.g., 
/i60Mm//i00Mm ~ 2 (see Table . For comparison with the 
sky confusion noise estimated from the ISO mission, we plot 
the HB90 formula to which our empirical relation is applied 
(see thick dotted line in Fig.jHJ. Although our formula solve 
the discrepancies of our estimations to some extent, there 
are still disagreements especially with the results for higher 
resolution missions. 

The HB90 formula was obtained from the analysis of 
the low resolution IRAS data at 100 iJ,m and assumed a 
constant power index for the cirrus power spectrum. In the 
case of the high resolution missions, since the sky confu- 
sion becomes sensitive to the local structure rather than the 
large scale structure, the calculation of the sky confusion 
strongly depends upon the power spectrum estimated for 
each patch and the power at the scale length correspond- 
ing to the resolution of the detector. Therefore, we should 
consider carefully the combination of the resolution and the 
power spectrum of the cirrus in the estimation of the sky 
confusion noise. In addition, the larger discrepancy in the 
bright regions for the ASTRO-F mission compared to the 
prediction from ISO observations can be explained by an in- 
crease in the spatial resolution, although the aperture sizes 
of two telescopes are similar (see the specifications of the two 
space missions in Table 0. We conclude that the sky con- 
fusion level predicted by the IRAS data from which HB90 
formula are derived is significantly overestimated in the case 
of the higher resolution missions. 

Generally the most important component superim- 
posed on the extragalactic background in the far-IR is 
the cirrus emission. However, at high spatial frequencies 
the Cosmic Far-IR Background (CFIRB) fluctuations 
may become dominant I jSchlegel, I 'inkbeiner fc Davis 1998| 
Guiderdoni et al. 1997 [Juvela, Ma ttila & Lcmkc 200d]L 
Therefore, in any estimation of the sky confusion noise 
using observational data in the dark fields should consider 
the fluctuation due to the CFIRB. By fitting the sky 
confusion noise over the mean sky brightness. Kiss et al. 
(2001) obtained CFIRB fluctuation of 7 ± 2 mjy at 90 
/xm and 15 ± 4 mJy at 170 /^m. After correcting for the 
contribution of the CFIRB in the estimation of the sky 
confusion noise, we obtained results similar with those of 
Kiss et al. Il2001t in the dark flelds (see the symbol in 
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Figure 8. Estimated sky confusion noise for the ISO mission. Upper and lower panels show the sky confusion noise at 90 fim and 170 
/im, respectively. The dotted line shows the sky confusion noise by HB90 (Helou & Beichman 1990). The symbols are the estimated 
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Figure 9. Estimated sky confusion noise for the ASTRO-F mission. Left and right panels show the sky confusion noise in the WIDE-S 
band (75 /jm) and WIDE-L band (140 ^tm), respectively. The symbols and lines are same as given in Fig. ISl 



circle with arrow in Fig. |H|at the mean brightness of ~ 1.5 
MJy/sr). Since the CFIRB fluctuations strongly depend 
upon the extragalactic source count model, we will discuss 
this issue in greater detail in our forthcoming paper [Jeong 
et al. 2004c l |Paper iT^ , in preparation]. 

4.2.3 Sky Confusion Noise for Various Separations 

Kiss et al. II200H analyzed the dependency of the sky confu- 
sion noise on other separations by the simple power expres- 
sion from ISO observational data: 

iV(qemin) = iV(e^in) X q< , (12) 

where q> 1 and 7 is a constant for a specific map. We ob- 
tained 7's for all patches and showed 7 as a function of mean 
brightness for each mission as given in Fig. 1131 As the sky 
becomes brighter, 7 becomes larger due to the prominent 



structure of the cirrus emission. Kiss et al. ll2UUlt obtained 
a much lower 7 in dark regions, but their values of 7 in 
other regions are similar to our results. This result can be 
explained by two possible effects: one is that the cirrus struc- 
ture observed by ISO is blurred by the instrumental noise in 
most of the dark regions and the other is that many extra- 
galactic point sources below the detection limit, i.e. CFIRB 
fluctuations, can remove the cirrus structure. If we only con- 
sider the component due to the cirrus in the dark fields, the 
values of 7 in the dark regions by Kiss et al. T^MTl are sim- 
ilar to our results. In most of the bright regions, the scatter 
of 7 shows the similar trend and this is probably caused 
by the relatively large difference in the spatial structure in 
each region. In the same mean brightness, 7's in SW band 
are larger than those in LW band because spatial structures 
should be prominent in SW band. In addition, since we use 
the simulated data, changing features of 7 in two wavelength 
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Figure 10. Estimated sky confusion noise for the Spitzer mission. Left and right panels show the sky confusion noise for the MIPS 70 
fim and 160 ^m bands, respectively. The symbols and lines are same as in Fig. |^ 
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Figure 11. Estimated sky confusion noise for the Herschel and SPICA missions. Left and right panels show the sky confusion noise at 
70 iim and 160 fira, respectively. The symbols and lines are same as in Fig. IHl 



have a similar shape. For the Herschel and SPICA missions, 
our estimations show that 7 slowly increases and the er- 
ror decreases compared with other missions, because of the 
much higher resolution than the other missions considered. 



4.2.4 Effect of Power Index a 

In this study, we assume that the structure of cirrus is in- 
dependent of wavelength. However, recent papers reported 
on enhanced dust emissivity at some medium-to-high den- 
sity clouds in LW band of Far-IR due to the presence of 
a cold dust component (T < 15K) (Cambrcs y et ai. 2001| 
|dei Burgo et ai. 2003| [Stepnik et ar'2003l . This result im- 
ply that the cirrus structure can be changed in LW band. 
Kiss et al. (I2UU3II suggested that the power index of the 
power spectrum also depends upon both wavelength and 
surface brightness due to the coexistence of dust compo- 
nents with various temperatures within the same field and 
cold extended emission features (usually, —2.0 < a < —4.0). 
Using the assumption that the sky confusion noise is propor- 



tional to the scale length (see equation 3J, we can estimate 
the sky confusion for different power indices. The ratio t/i of 
the sky confusion noise with the power index of a-f e to that 
with the power index of a can be defined as 



N{a + e) 
N{a) ' 



(13) 



where e is the contribution to the power index from any 
other structure in the power spectrum. In this calculation, 
we fix the power at the scale length of the resolution limit of 
the map (~ 6.1 arcmin) and wavelength at 100 from the 
assumption that the power over this scale is not affected by 
the extra components proposed by Kiss et al. (12003^ . Table 
2]hsts the ratio of the sky confusion noise for the different 
power indices for each space mission covering power indices 
of the power spectrum on the cirrus emission. Since the fiuc- 
tuation at smaller scales is more sensitive to the power in- 
dex, the sky confusion noise is much more dependent upon 
the power index for the space missions with higher resolu- 
tions. As seen in Table|3 our estimated power indices in the 
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bright regions (a > 3.3) are somewhat higher than those in 
low density regions {a < 2.8). From the recent Spitzer ob- 
servation, Ingalls et al. 120U4t obtained the power index of 
-3.5 at 70 j-im in the Gum Nebula. Therefore, if this varying 



Table 4. Ratio ip of the sky confusion noise for the different 
power indices. 
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power index is not so large, it does not affect severely the 
final sensitivity values. 



5 PHOTOMETRIC MEASUREMENTS OF SKY 
CONFUSION NOISE 

In Section 2] we estimated the sky confusion noise by the 
fluctuation analysis. The sky confusion noise should affect 
the source detection efficiency, causing a deterioration in the 
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detection limit. In this section, we obtain the measured sky 
confusion noise by carrying out photometry on realistically 
simulated data. 



5.1 Source Distribution 

The distribution of sources per unit area on the sky can 
be described as a function of the flux density and depends 
upon both the spatial distribution of the sources and their 
luminosity function. For simplicity, we assume the number 
of sources whose flux is greater than flux F, n(> is a 
power- law function of F, 



n(> F) = no(> Fq) 



(14) 



for F„, 



F < -Fmax, where no and Fo are normalization 



constants for number of sources and for flux, respectively, 
-Fmin is the minimum flux, Fmax is the maximum flux in the 
source distribution. 

The source confusion caused by the overlap- 
ping of adjacent sources mainly depends upon the 
source distribution and the beam proflle (|Condon 19741 
IFranceschini et al. 1989|l . Source confusion becomes impor- 
tant as the observation sensitivity increases since there are 
usually more faint sources than brighter ones. Currently 
favorable source count models require strong evolution in 
order to fit the ISO data from mid- to far-IR, the SCUBA 
data at sub-mm wavelengths, and the Cosmic Infrared 
Background (CIRB) at 170 fim llOliver et al. 19971 

IKawara et al. 19981 



Smafl, Ivison fc Blain i997| 

Hughes et al. 199 8| lAussel et al. 19991 |Puget et al. 1999 1 
'Efstathiou et al. 2000l [Serjeant et al. 20001 

IMatsuhara et al. 200ol 

In our study, we use a simple source 



Lagache et al. 2000| 
Scott et al. 2nn'2r~^ 



distribution for the purpose of investigating only the effect 
of the sky confusion. We will discuss the source confusion 
with more realistic source count models in the forthcoming 
paper. In order to avoid the contributions from any source 
confusion itself, we assume rather sparse distribution of 
sources. However, the estimate of detection limit becomes 
rather uncertain, if there are too few sources. Therefore, we 
have employed a model for the n{F) utilizing a distribution 
with two slopes, uj — 1.0 for bright flux region and u — 
0.3 for faint flux region (see Fig. I14L in order to derive an 
accurate value for the sky confusion limits without source 
confusion effect. Since the sky confusion noises in the SW 
bands are much lower than those in the LW bands, we set 
different normalization constants and minimum flux values 
Fmin, i.e., Fmin = 0.001 mjy and no(> Fq) = 3 in the SW 
band, Smin ~ 0.1 mJy and no(> Fq) = 10 in the LW band, 
where _Fo is set to be 100 mJy (see Fig. I14II . 

5.2 Source Detection 

We generate images including point sources convolved with 
the beam profile of each mission using the source distribu- 
tion described in Section 15.11 Fig. 1151 shows the simulated 
images for the various missions considered. As the detector 
pixel and the beam profile become smaller, more sources and 
smaller structure in the cirrus emission appear. 

We carried out aperture photometry on the sim- 
ulated images using the SExtractor software t;2.2.2 
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Figure 14. Source distribution in the SW band and LW band. 
We use different slopes [ui = 1.0 and ui = 0.3) for the power law 
source distribution at the boundary flux of 10 mJy in order to 
reduce the effect of the source confusion. 



IBertin fc Arnouts 199611 . There are several parameters to 
be fixed to perform the photometry, but the most influen- 
tial parameters are the size of a background mesh for esti- 
mating background level and the threshold for the source 
detection in this aperture photometry. In order to optimise 
for better reliability of the detected sources and reducing 
the rate of false detection, we make trials by changing two 
parameters. Finally, we set the size of the background mesh 
to be 2.5 times of the measuring aperture, and the detection 
threshold as 4a. The final detection limit is determined by 
the minimum flux of detected point sources. We found that 
the detection limits determined from 4a criteria are consis- 
tent with the 4 times of sky confusion noise measured from 
the fluctuation analysis. Note that our sky confusion noise 
estimated from the fluctuation analysis is a la fluctuation. 
In Fig. 1161 we compare the detection limit by photometry 
with the sky confusion noise for each mission. For the ISO 
and ASTRO-F missions, the results from photometry give 
relatively higher detection limits than the theoretical esti- 
mations via fluctuation analysis. This trend results from the 
larger detector pixel size compared to the FWHM of the 
beam proflle. The large detector pixel size of the ISO mis- 
sion significantly degraded the performance of the detection 
of the point sources (e.g., the left panels in Fig. ll6ll . 



6 SUMMARY AND DISCUSSION 

Based on the observed 100 /im dust map and the models of 
a dust spectrum, we generated high resolution background 
maps at wavelengths ranging from 50 to 200 /xm. Using these 
simulated cirrus maps, we estimated the sky confusion noise 
for various IR space missions such as ISO, Spitzer, ASTRO- 
F, Herschel and SPICA. Since we have the observational 
results only from ISO, we compared the results of our simu- 
lation with the ISO data. We found that the sky confusion 
noise estimated with our simulated maps are consistent with 
the ISO results. However, in the dark fields the sky confusion 
noise is more weakly dependent upon the beam separation 
parameter than in the bright fields in the case of the ISO 
observation. We conclude that this is due to the fact that 
the instrumental noise dominates in the dark regions or al- 
ternatively, the CFIRB fluctuation is more important. We 
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Figure 15. Simulated images including point sources in the LW band for ISO (upper-left), ASTRO-F (upper-right), Spitzer (lower-left), 
Herschel and SPICA (lower-right) missions. The mean brightness of the cirrus background is 2 MJy sr~^ at 160 fira. 



also found that the sky confusion predicted from the IRAS 
data is significantly overestimated in the case of the large 
aperture telescopes, except for the dark fields. 

We have confirmed our results through a realistic sim- 
ulation. We performed photometry on simulated images in- 
cluding point sources with a sparse source distribution in 
order to avoid the effects of confusion due to crowded point 
sources. The detection limits obtained from the photometric 
analysis agree with the sky confusion noise estimated using 
fluctuation analysis except for ISO and ASTRO-F. The dis- 
crepancies for these missions are due to the large detector 
pixel size compared to the FWHM of the beam size. 

The mean brightness of the cirrus emission 
usually decreases with increasing Galactic latitude 
l |Boulanger fc Perault 1988| . In order to estimate the 
detection limits as a function of Galactic latitude, we 
derived a simple formula for each wavelength band. Because 
the cirrus emission is extremely strong near the Galactic 
centre, we excluded the Galactic latitudes |fe| < 10°. Fig. 
1171 shows the detection limits as a function of Galactic 
latitude. The detection limits for all missions appear to 
saturate beyond b ~ 30° . 

Fig. 1181 summarises the final detection limits for point 
sources at mean and low sky brightness regions due to the 
Galactic cirrus. In addition, we also plot the currently es- 
timated 5a detection limits for sources of each mission. 
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Figure 17. Detection limits due to the Galactic cirrus as a func- 
tion of Galactic latitude. The two line plotted for each mission 
are for the SW band (lower line) and the LW band (upper line). 



The detection limits only take into account the instrumen- 
tal noise. The instrumental noise for ASTRO-F mission is 
explained in detafl in Jeong et al. (2003; 2004a; 2004b). 
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Figure 16. Estimated detection limit by photometry. Figures show the detection hmit and 4 times sky confusion noise estimated from 
the fluctuation analysis for the ISO and ASTRO-F missions (left) and Spitzer, Herschel and SPICA missions (right). Upper and lower 
panels show the results for the SW band and LW band, respectively. 



The integration time is 500 sec for the Spitzer mission 
{Spitzer Observer's Manual^) and 1 hour for the Herschel 
mission iPilbratt 200311 . As shown in Fig. 1181 sky confusion 
almost approaches the detection limit in the LW band of the 
ASTRO-F and Spitzer missions. Although the sky confusion 
does not severely affect the detection limits of Herschel mis- 
sion, it can affect the detection limit of the SPICA because 
it will have a large aperture telescope cooled to very low 
temperatures in order to achieve exceptional sensitivity in 
the far-IR (see Nakagawa 2004 for the detailed information 
of the SPICA mission). 
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